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Abstract 

We have constructed a one dimensional exactly solvable model, which is based on the 
t-J model of strongly correlated electrons, but which has additional quantum group sym- 
metry, ensuring the degeneration of states. We use Bethe Ansatz technique to investigate 
this model. The thermodynamic limit of the model is considered and equations for different 
density functions written down. These equations demonstrate that the additional colour 
degrees of freedom of the model behave as in a gauge theory, namely an arbitrary distribu- 
tion of colour indices over particles leave invariant the energy of the ground state and the 
excitations. The 5-matrix of the model is shown to be the product of the ordinary t — J 
model 5-matrix and the unity matrix in the colour space. 
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1 Introduction 



Since the discovery of high-Tc cuprate superconductivity the one dimensional physics of 
strongly correlated electrons has been in focus in many publications |l|. The Hubbard Q 
and t — J models Q are such examples , motivated in part by the fact that high-Tc com- 
pounds display antiferromagnetism in the absence of doping. The t — J model was proposed 
by Zhang and Rice and describes strongly correlated electrons with anntiferromagnetic 
exchange interaction. 

The interest in one dimensional physics grew after Andersons claim |Q] that two-dimensional 
systems may have features in common with one dimensional systems. In addition it should 
be mentioned that powerful methods in ID such as bosonization, 2D conformal field theory 
and in particular the Bethe ansatz technique allow the detailed study of such systems. 

The t — J model may be used as well for heavy fermion system ^. 

At the supersymmetric point J = 2t the t — J model becomes exactly integrable 
m, |l^, 11, 12 1 because the Hamiltonian can be represented as a graded permutation in a 
super algebra of two fermions and one boson. 



In [03, 14, 15| we developed the technique for construction of a family of spin chain 



Hamiltonians and their fermionic representations, which have the same energy levels as 
some basic model{XXZ, Hubbard, t — J or others) but with huge degeneracy as a result 
of an affine quantum symmetry added to the basic model. We called this procedure an 
affinization of the model. 

The first example of this type of model was constructed in ||l^ , giving rise to the Hubbard 
Hamiltonian in the infinite repulsion limit. 

In 1 14 1 we have fermionized the simplest examples of this newly defined family of models 
and have shown that it leads to extensions of one-band Hubbard Hamiltonians. The re- 
pairing mechanism introduced by Yang |l9| was found in one of examples in addition to 
other exactly solvable Hubbard models with superconducting ground state |2^, |2^, 
24|. The essential property of this extensions is the fact that, besides the ordinary electron 
hopping and Hubbard interaction terms, they contain also bond-charge interaction, pair- 
hopping and nearest-neighbour interaction terms. In |l^] the SU{N) affinization of t — J 
model was carried out, giving rise a model where the spin-spin coupling term consists 
of interaction between the total spins(i.e. the sum of the spins of all band) at nearest- 
neighbour sites. The presence of the affine symmetry, which ensures the degeneracy of 
levels exponentially proportional to the length (area) of the space, might lead to a new type 
of string theory. 

In this article we define an extension of t — J model such that an affine quantum group 
symmetry is present, and we use the Bethe ansatz technique to solve the model. We 
find the S-matrix excitations on empty background, the ground state and construct the 
thermodynamic limit of the model. As one might expect, the S-matrix of the excitations on 
empty background consists of the S-matrix of the ordinary t — J model multiplied by the 
unity matrix in the additional space of "colours". Therefore the Bethe equations are not 
different from ones for t-J model, but the rapidities presented in equations correspond to 
particles with the arbitrary colours. The degeneracy of the corresponding n-particle states 
come from arbitrary partitions of the colour indices over particles. The same is true for 
the ground state. The situation is exactly as in gauge theories if we distinguish the states 
which differs by pure gauge transformations. All this results are presented in the Section 5 
of the article. 

The thermodynamic limit of the model with corresponding equations are represented 
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in the Section 6, where we also shown that the S-matrix of our model in an arbitrary 
background is equal to ordinary t — J model S-matrix multiplied by the Kronecker symbols 
over the additional colour indices. 



2 Quantum group invariant Hamiltonians for reducible rep- 
resentations 



Let V = ®iLiVx- be a direct sum of finite dimensional irreducible representations V 



A,: 



of quantum group UqQ |32, We denote by V(xi, . . . ,xn) the representation with 



spectral parameters Xi of the corresponding affine quantum group Uqg |34] 



M 



V(xi,...,x^) = 0NA,0VA,(xi), (1) 

i=l 

where all the \'\.{xi) are M nonequivalent irreps and N;^, ~ C^' have dimensions equal to 
the multiplicity of \\^{xi) in V(xi, . . . , xtv)- Note that Y^fLi = The ' over the tensor 
product signifies that UgQ does not act on N;^-®V;^. (xj) by means of co- multiplication A 
but instead acts as id (8" g. 

In [^] the general matrix form of the intertwining operator 

H{xi, . . .,xn) ■■ 

V(xi, . . . , xn) (g) V(xi, . . . , xn) V(xi, . . . , xn) (g) V(xi, . . . , xat), (2) 
[H{xi,...,XN),A{a)]=0, yaeUqg 

had been written using the projection operators 

XS = \a){b\ (3) 

Here the vectors \a) span the space V. In accordance with the decomposition (||) we will use 
the double index a = (n^, aj), i = 1, . . . ,M where the first index rij = 1, . . . , A'j characterises 
the multiplicity of \^x. and the second one a, = 1, . . . ,dimV;^. is the vector index of V;^. . 
Then the intertwining operator (0) is 



M 



i,j=l \ni,nj,mi,mj ai,aj 



where the i?-matrix 

Rvx^(S)VxX^i/xj)\ai) ® \aj) = Rij''j''j{xi/xj)\a'i) \a'j) 



a' ,a' 



is the intertwining operator between two actions of affine quantum group Ugg on V^^ ® , 
which are induced correspondingly by co- multiplication A and opposite co- multiplication 



A [mH, 
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Aij and Bij, Bu = in are arbitrary matrices. In general, H{A,B) depends on de- 
formation parameter q of quantum group, which is included in the i?-matrix. Note that 
-?^Va .iX)Va . (^j/^i) does not depend on q and is identity only if Aj or Xj are trivial one- 
dimensional representations. So, in the special case, when the only nontrivial i?-matrixes 
in (^) are between representations, one of which is trivial representation, the expression of 
H{A, B) doesn't depend on q. Then H(A, B) commutes with the quantum group action for 
all values of deformation parameter. In the following we consider only this case. 

Following |l^ we can from the operator H construct the following Hamiltonian 
acting on W = 

L-l 

n=Y. Hu+i, (5) 

i=l 

where the indices i and i+1 denote the sites where H acts non-trivially. By the construction, 
7i is quantum group invariant: 

[H,A^-i(5)] =0, V5Gf/g5 

Let us define the projection operators on V for each class of equivalent irreps (Aj,Xi), 
i = 1,...,M 

M 



i=l 

Their action on W is given by 



I 

1=1 

It is easy to see that these projections commute with Hamiltonian Tl and quantum group 

[Q\n] = o, [Q\US = (6) 

Denote by Wp^...p^ the subspace of W with values pi of on it. Then we have the 
decomposition 

W= Wp,..p,, (7) 

P1,...,PM 

pi+...+PM=L 

Let V'^ be the linear space, spanned by the highest weight vectors in V: 

where Vy € \x. is a highest weig ht vector. We also define w° := V° The space W" 
is 'K-invariant. For general q the action of Uqg on W'^ generate whole space W. Indeed, 
the Uqg-action on each state of type v^. . . v^^ generates the space V;^.^ . . .<0 ^^iL ' 
because the tensor product of finite dimensional irreducible representations of an affine 
quantum group is irreducible |35|. 

^Here and in the following we omit the dependence on Xi 
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Consider now the subspace = W*^ n '^pi...pM- According to (0) we have the 

decomposition 

w°= e w°^,„p^^. (8) 



PI, ■■■,PM 

pi+...+PM=L 



Note that 



L 

Pi--- PM 



Let us define by Hq the restriction of H on Wq: TCq := H|wo- It follows from (^) that 
Hamiltonians Tl and T^o have block diagonal form with respect to the decompositions 
and (|8|), respectively. Every eigenvector w^^^ G ^pi-.-pm '^ith energy value E^p^ ^^^^ 
gives rise to an irreducible C/q^-multiplet W^^^ of dimension 



M 

dimW„^^...^^^ = n(dimV,JP^ (9) 

k=l 



On W„ the Hamiltonian Ti is diagonal with eigenvalue En ■ In particular, in the 
case when all V^^ are equivalent, the degeneracy levels are the same for all E^^^ and 
are equal to (dimVA)^. 

Now, let us assume we know the energy spectrum Eap^^^^^ for TCq. Then the statistical 
sum is given by 

ZHoiP)= E E ^M-PEap^,,,pJ, (10) 

Pl + ...+PM=i 

and it follows that the statistical sum of 7i has the following form: 

M '^Pl---PM 

Zn{P)= E n(dimV,Jf'= E ^M-PEn,,...,J- (H) 

Pl,---.PM k=l api...pj^j=l 
pi+...+PM=L 

So, if the underlying Hamiltonian Ho is integrable and its eigenvectors and eigenvalues can 
be found, then we know these for Ti. too. Acting with the quantum group on all eigenvectors 
of an energy level of Hq one obtains the whole eigenspace of 7i for this level. 



Multi-band t 

J = 



J model with vanishing spin-spin coupling 



Let us consider here the quantum group Ugsl2- We choose V = Vq © Vj for decomposi- 
tion (|l|), i.e. we take a direct sum of the trivial spin-0 and the 2j + 1-dimensional spin-j 
representation of Uqsl2- The ii-matrix in the second term in (^) does not depend on q and 
spectral parameters Xi and coincides with the identity, as it was mentioned above. So, using 
(^ and (^), we obtain the following Hamiltonian 



L-l 

nit,Vi,V2) = J2 
1=1 



2j+l 

E i-^iO-^i+^p + A'i+loA'ip) + FiXjQXj+io 

p=l 

2j+l 

+V2 E ^ip^i 
p,p'=l 



(12) 
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The Hamiltonian Ti = J^i^a+i constructed from the operator H = Ha^i, where H 
can be written in the matrix form 



H 



1 ^1 








\ 








-t ■ id 








-t ■ id 








V 








V2-id J 



(13) 



The projection on the highest weight space coincides with the constructing block of the 
XXZ Hamiltonian in an external magnetic field. This implies that the restriction of ( |l2|) 
to the space W'' is 



no{t,Wi,W2 



nxxz{t,A,B) 



L-1 



i=l 



B 



where 



A 



2t ' 



B = -{Vi-V2) 



(14) 



(15) 



For the special case Vi + V2 = 7io gives rise to the free fermionic (or equivalently XY) 
Hamiltonian (A = 0). 

The projection operators X^ are expressed through the fermionic creation- annihilation 
operators as follows 



X-O — r- V 



v.P - Vr^ 

XiP = ni^pV = Vni^p, Xi^o = i^-ni)V = V{l-ni) (16) 
Here we introduced the projection operator which forbids double occupation on all sites 



'p=YiVi, = n (1 - ^i,p^i,p') 

p^p' 



i=l 



and the total particle number = J2p''^i,p site i. 

After the substitution of the fermionic representation (|l^) into ([l^) we obtain 



L-l 

n{t,Vi,V2) = rY. 

1=1 



2j+l 

* E (^tp^i+^,P + ^t+l,p^i,p) + 
p=l 

VniUi+i - Vi{ni + rij+i) + Vi] V, 



(17) 



where V = Vi + V2. The chemical potential term —yiJ2i=ii^i + "-j+i) commutes with 
Ti. and can be omitted. So, up to unessential boundary and constant terms ( |T7| ) is a 
multicomponent t — J model with vanishing spin-spin coupling (J = 0) 



L-l 



i=l 



2j+l 

E i4,pCi+i,p + cii,pCi,p) + Vuim+i 
p=i 



L 2j+l 
+ E E 

,p'''^i,p'^i,p' 7 

(18) 



p,p'=l 



where the infinite Hubbard interaction amplitude Up^p' = +00 between p and p' bands 
excludes sites with double and more occupations. It follows from the above considerations 
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that this model has energy levels which coincide with the levels of XX Z Heisenberg model, 
but that the degeneracy of the levels is different. 

For vanishing density-density interaction V = Q the Hamiltonian ([l^) describes the 
infinite repulsion limit of the multi-band Hubbard model. Thus, according to (ITsI ) A = 
and it has the energy levels of the free fermionic model. 



4 Multi-band extension of the t — J model with afRne quan- 
tum group symmetry 

In this section we consider Hamiltonians which have the same energy levels as t — J model 
but have affine quantum group symmetry. Because each site in ordinary t — J model has 
three states, one should for this purpose take direct sum of three spaces. Let 



V = Vo e Vj- e Vj 

Recall that the t — J model is given by 

X! (Cj>Cj+i,<T + cf_^i^^Ci^a) + JSiSi+i + Vriim+i 



(19) 



L-l 

nt-j{t,j,v)=vJ2 
1=1 



a=±4 



V, (20) 



where c+, are creation-annihilation operators of spin-i fermion, S = X^o-o-' C^^aa'Ca' is 
the fermionic spin operator and V = ni^i(l ~ ''^i,T"'«,J.) forbids double occupation of sites. 
We rewrite it in terms of Hubbard operators X^, where a,b = 0, 



L-l 



1=1 



+ Yl ia<T'j + v)xcXi+c', 



a,a'=±\ 



(21) 



Let us now look at the general expression (^ of intertwining operators Hij acting on the 
space ([l9|). For convenience we make index change in the following way. The two spin-j 
representations we use are denoted by o" = zti. The intrinsic index in each V"/'^^ is denoted 
by A;, /c = 1, . . . , 2j + 1. So, instead of (n^, Oj) in (^) we have (a, /c), if i corresponds to spin-j 
multiplet. Because the spin-0 singlet is one dimensional and single, we just use for it the 
index 0. The non-equivalent irreps in (|9|) are V^"^ and Vq and, as mentioned above, the 
i?-matrix for two such representations is the identity. After performing the first sum in ^ 
over non-equivalent multiplets we obtain 



H{A,aMM) 




criCTj Z / (o'2,k) 

k,k' 



X, 



+ 



bi-X. 



{a,k) 



^{a,k) + ^2 • X(a,k) 



X, 



(a,k) 



(22) 
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To implement the restriction Hq{A, a,bi,b2) of this operator on the highest weight space 



one just should eliminate the sum over k, k' and put k = k' = 0. Comparing (22) and 
it follows that the expressions coincide if one chooses 



a = bi = b2 = -t A-Z = J/2 



A" 



[aa') -J + V 



and choose the other values of A'^aXZl equal zero. 

So, the Hamiltonian 7Y(A, a, 6i, 62) corresponding to (22) with these values of parameters 
gives rise to a t — J model ( pO|) on the highest weight space. According to the previous con- 
siderations it will have the same energy levels as t — J model, but with different degeneracy. 
Recall that for J = 2t the t — J model is "supersymmetric" and integrable. 

We express the Hubbard operators in terms of multi-band fermionic creation-annihila- 
tion operators as follows 



Vc: 



Y 



(23) 



k+ k -p _ T)(M+ k jr{^,k) 



«(-o-,fe) 

Here as before we used the projection operator, which forbids double occupation on all sites 

L 



i=l 



(cr,fc)7^(cr',fc') 



Now, we can write down the Hamiltonian (|5|) in terms of multi-band fermions, substituting 
(p^) into (22). We obtain in this way the multi-band generalization of ^ 



L-l 

nit,j,v) = vJ2 

i=l 



2j+l 



k=l a=±\ 



V, (24) 



Here k is the band index, and S = J^k^'^^ ^ — Ylik"^^ total spin and total particle 
number operators. It's easy to see that we have conservation of the particle number operators 
Y^^a- — "^^ foi' the all colours k. 



5 Bethe ansatz for the t— J model with affine quantum group 
symmetry 

The goal of this section is to apply the Bethe Ansatz technique to our model and derive the 
corresponding Bethe equations for the excitations. 

After making some trivial Pauli matrix calculations one can represent Hamiltonian ( p4|) 

as 



L-l 

n{t,j,v) = vY. 

i=l 



2j+l 



E cf+cf,,4+ ^cf+i,, + (l^ - J)nini+i 

k,k' ,a^T 



V. 



(25) 
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Due to the conservation of the particle number operator n'^, and according to the coordinate 
Bethe Ansatz we look for eigenvectors of (^), corresponding to N fermions of 2j + 1 bands 
in the following form 

1^ >= E E - E i^''---'Hxiai,---,XNaN)4\t,--4%tjo >, (26) 

where |0 > is empty vacuum state. 

The eigenvalue equation >= E\^ > in the one particle sector 

- t{^''{x + ^''{x, a)) = E^''{x, a) (27) 

gives us 

E = -2tcosp (28) 

after substituting of the plane wave function with momentum p into (|27| ) . 

The eigenvalue equations in the two particle sector allows us to fix the energy of the 
state as a sum of two one particle energies, as well as the two two particle scattering matrix 

a' a' k' k' u h 

^a\cT2' kik^iP^-'P'i)- choose the antisymmetric wave function '^{xiai^X2(T2) as 

V''i''=2(xiai,X2CT2) = A^i'=2(piai,p2fT2)e*^Pi"i+P^"2) - ^'=2'=HP2<Ti,pia2)e^(f2"i+fi"2) (29) 
for X2 < xi and 

V''=i''=2(xi^7i,X2a2) = A^2'=HP2^T2,PitTi)e*^Pi"i+*'^"^) - ^'^^'^^ (rcts, P2tTi)e'^^'2"^+^^"^^ (30) 
for xi < X2- 

The continuity condition at xi = X2 should be imposed: 



Use of the Hamiltonian ( [24| ) and the most general form ( [2q ) of the eigenfunctions ^, the 
eigenvalue equations can be written as 



t 



i^^^^^{xi + lfTl,X2fT2)(l - 5^^ + l^^^)+i)^^^''{xi - lfTl,X2fT2)(l - 5xi,X2 + l) + 
'(xifTi,X2 + lo-2)(l - 5xi^x2-l) + t + ^^^^^{xiCFi,X2 - lcr2)(l - ^aii-l.xa) + 

'iJ5\xi-x2\,i^^'^^{xia2,X2ai) + +{V - J)5\^^^^^iii)^'^^{xiai,X2a2) 



= Ei,^'^^{xiai,X2a2) (32) 

The terms (1 — 5) appeared near the hopping terms as a result of projective operator 
V, preventing double occupancy of the sites. 
Two different cases can be considered: 



I. |xi — X2I > 1. In this case eigenvalue equation ( |32| ) is reduced to 

-t[ilj''''"'{xi + lai,X2a2) + i^^^^^{xi - lai,X2CJ2) + i^^^^^{xiai,X2 + 1^x2) + 

+ij^'^^{xiai,X2 - la2)] = E^^^^^{xiai,X2CJ2). (33) 



After substitution of the expression (p9|)-(|30|) for the plane waves into (|33| ) and some 
simple calculations, the spectrum of two particle state can be found. 

E = -2i(cospi + cosp2)- (34) 



9 



II. I^i — X2I = 1. Without loos of generality one can take X2 = xi + 1. Then the 
eigenvalue equation reduces to 

-t[i;''^''^{xi - lai,X2a2) + i^'''''%xiai,X2 + la2) + 2,711;''^''^ {xia2,X2C7i) + 

+ {V - J)^^'^^{xiai,X2(T2) = Eij^'^^{xiCJi,X2(J2) (35) 



By substitution of (^) and use of continuity conditions ( pl| ) one can express the ampli- 
tude A^'^ {p2 ai , (T2 ) after scattering through A'^^ (pi cri , p2 ^2 ) and A^^^^{pia2,P20'i) 
before, and, therefore get R-matrix of the model: 

A''''HP2ai,pia2) = R''''"^'lf'A'^>'^{pia[,p2a2) (36) 

where i?-matrix is the product of ordinary t — J model i?-matrix multiplied the per- 
mutation operator in the /c-index space. 

(Ticr2,A;jA;2 ^ ' ^ ' 

At the supersymmetric point 2 J = t, V = — J/4 



Ai - A2 + z 



In (|38|) A = i cot I is the rapidity, I = SZl^al ,and P = (5^1 (5a 



The scattering matrix SJ^J;^^^-^^^ will be defined multiplying the i?-matrix by permuta- 
tion operator P in the spin((T) and colour(/c) spaces. 

The exact integrability of the model is connected with the fact that the S'-matrix should 
fulfill the Yang-Baxter triangular relations 

5:ii£6(A, - A2) • 5:^3^3 _ . 5g|f3(A2 - A3) = 
^iia(A2 - A3) • ^x}jSl^' - • 4Z^^'' - ^''^ 

and constraints on t, J, V are imposed just by these equations. 

Consider now N itinerant electrons in a box of L sites with periodic boundary conditions. 
If we successively make change of position of an electron and its neighbouring electron in 
a chain, each interchange produces a scattering matrix and when the particle comes back 
to its starting position from the other side, we will have the cyclic product of S'-matrices, 
which is called the Transfer Matrix: 

^i('^i) = 'S'i,i+i(Aj — Aj+i)....5'j,Ar(Aj — Aat) • Sj^i{Xj — Xi)....Sjj^i{Xj — Xj~i) (40) 

Here we skip the matrix indexes a and k while the hat on S means the operator in that 

space. The periodicity means that the Transfer Matrix has to be diagonal for all j = 1 

with the eigenvalues exp{ipjN), or, on terms of rapidity Xj 

^ip,N ^ ( --J ' y - ] (41) 
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The matrix ([40| ) is the trace of so cahed monodromy matrix, which by definition is the 
product of the S-matrices without taking trace. Hence, the monodromy matrix can be 
considered as a (2 x 2) (g) {{2j + 1) x (2j + 1)) matrix in the spin a and colour k spaces. 
Let's remember now that this operator is a unity operator in the colour space. We will 
not describe here the details of the algebraic Bethe Ansatz(see e.g. [^0|), but already now 
it is clear, that because the monodromy matrix of our model is the monodromy matrix of 
the ordinary t — J model multiplied by the unity matrix in the /c-space, the generalization 
of the algebraic Bethe Ansatz to the present model give rise to the same equations as the 
ordinary t — J model. 

Specifically we get: 

where L is the number of lattice sites, A^-is the number of electrons and AI is the number 
of spin down electrons. 

We see that colour is disappeared from the equations, which means that we can make 
arbitrary partition of colour charges on the state of A'^ particles in a chain and all wave 
vectors will become eigenvalues of the Hamiltonian, provided that their A's fulfill the Bethe 
equations. 

Equations (p2[)-(^3|) are Lai's Q form of Bethe equations, written on a basis of empty 
background. 

The total energy is given by 

^ 1/2 

However for the construction of the thermodynamic limit of our model it is more con- 
venient to use the Sutherland's form of the Bethe equations, which is equivalent to Lai's 




equations |11, In Sutherland's representation one start with a ferromagnetic pseudo- 
vacuum state with all spins up and consider excitations as A^/^-holes (empty sites) and M 
spin down electrons. 

In Sutherland's representation we have 

n ^tei^. "-'-■-^ (^^) 

The equations ( p5| ) have real and complex solutions. The complex solutions are in a form 
known as strings, which may be found by fixing and M (they are conserved quantities) 
and letting the lattice size L ^ oo. Following 0, |6| one finds complex solutions in the 
form 

Xj = \ + i{n -|- 1 — 2j), j = 1, ...n 

A^ = \ + i{n-2T) r = l,...n-l, (46) 
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for arbitrary n. In a finite box these string solutions are not exact but as in ||8|, |2^, we will 
assume that the corrections are small in power of L and hence vanish in the thermodynamic 
limit. 

We would like to write the equation for the centrum of the strings of length n, A^, (which 
are the real variables), let their number is M and rapidities of spin down electrons. 

Substituting ( [45D into Lai's form of Bethe equations and taking the logarithms from the 
left and right hand sides, one can get 



yn oo Nm M \n _ A 

m=l 7 13=1 



and 



oo Nn A , _ \n 

E E — —) = ^^Jp (48) 



1 1 n 



where is the number of strings of length m, 9{x) = 2tan ^{x) and 

a> = S ^(^) + 2^(FI=fe) + - + 2^(^) ri^m 

""^"^ \ 20(f) + 20(f) + ... + 20(^^) if n = m ^ 

and Jj3 are integers(half integers) appearing after choice of branch of the logarithm. 
The integers M, and M satisfy the relation 

oo 

M' = M + E (50) 

m=l 

Any solution is defined by a set of I" and Jp. All possible integer (half integer) values defines 
the states called vacancies. The vacancies may be occupied by particles of colour k or not 
occupied at all. In a case of occupied vacancies one can mark the corresponding and Jg 
by colour index k as for particles, I^'^ and J^. When the state is empty, the corresponding 
integer will be mentioned as and Jp. 

One may calculate the number of string states (^) taking account also colour degrees 
of freedom and will get {2j + 1)'^. 



6 The thermodynamic Umit of the Bethe Ansatz equations 

In the thermodynamic limit L — > oo the solution becomes densely packed (the difference 
of two neighbour solutions is of order 1/L, Xj+i — A" = 0{1/L)) and one can introduce 
density functions of the states and pass from sum's to integrals in the equations 
Let's define the particle and hole densities as follows: 



p^'^iX) = hm — — pl{\) = hm 



»00 L(X ri,k — Aj-n.fc) L— »00 L(^\jn — ) 
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and the sum of the hole and particle densities defines the density functions of vacancies. 

2j+l 2j+l 

P?=T.Pp'' + Pt at=Y.^P+^h- (52) 

k=l k=l 

The passage from sum to integral is straightforward and in the thermodynamic limit 
the Bethe equations (|47|)-(H8D become integral equations: 



pUX) = /n(A) -J2 dX An,miX " A )p^(A ) - / dAfniX - A)cTp(A) (53) 



and 



roc 

at{A) = Y. dA/„(A-A)p^(A) (54) 



fn{X) and An^m{X) in equations (p^)-(|53) are defined by 



and 



/„W = f^f>=i-^ (55) 

ZTT dX vr + 



A ^ j /|n-»n|(^) + 2/|n-m|+2(A) + ... + 2/„+m_2(A), if Tl^m , . 

"'"^ I 2/2(A) + 2/4(A) + ... + 2/2„_2(A), if n = m ^^^^ 



Following Yang and Yang |^^, Takahashi [^] and Gaudin [^^, we can write down 
the thermodynamic equilibrium equations. For the ordinary t — J model it was done in 

The conserved quantities of our model are: 

• the energy (the expression in ferromagnetic background slightly differs from the Lai's 
expression (p^) 

/ oo 2j+l \ 

E = Lil-47rY.Yl J dXfn{X)p^/{X) , (57) 



oo 2j+l 

EE 

n=l k=l 

• the number of different particles 



N'' = L I 



j dxf2^p-/{X) + jdKaliK)^ (58) 



and the magnetisation 



= t ( 1 - E E (2^ - 1) / ^Vp''(A) - E / dkalik) \ . (59) 



oo 2j+l 2j+l 



n=l k=l k=l 



Correspondingly, for the free energy F one can write 



2j+l 

F = E-Y, PkN^" - BS, - TS, (60) 

k=l 
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where sue the chemical potentials of particles, B is the external magnetic field and T is 
the temperature. We will calculate the entropy 5 in a standard way, as the logarithm of 
the number of possible states of the model in the interval dX. 



N{X, dX) = e^(^)"'^ 
By Stirling's formula one finds from (|6l| ) 



U'^^iiLpkdxy.{LphdXy. 



(61) 



S = L 



171=1 



2j+l 



+ / dA 



k=l 
2j+l 

at log at - ^ ap log cjp - an log an 

k=l 



(62) 



To obtain the equilibrium equations one should minimise the free energy over density func- 
tions of the particles and holes. Puting 5F = for variations of the density functions 
Pp''', p^, ap,ah one gets after some algebraic transformations an infinite set of integral equa- 
tions for the densities. If one defines the excited energies as usual 



n n-,k 
Pt - Pp 



a" 



o-t - o-« 



e T 



and 



we obtain 



n,k 



-47r/„(A)-(2n-l)S + /ifc 



(63) 
(64) 



e-'f^iX) = e"^'" + I dA'A„,„(A - A') log (1 + 



m=l 



-T J dA/„(A-A)log(l + e— 

for n = 1, ....oo, and 

ei = -pk-B + TY^ J dA/„(A-A)log(l + e— 



(65) 



(66) 



n=l' 



It is seen from these equations for densities and excitation energies that if all particles 
have the same chemical potential p, = p^, the solutions will be independent of k. Corre- 
spondingly, for minimum value of the free energy one can has 



- = {l + B- p)-TY, J dX fmiX) log (l + e^ , 

m=l 

= (1 - 2p) - eA(A = 0) 



(67) 



in correspondence with [^, 26 1 

In order to analyse the ground state further one should take the limit T ^ in equations 



3|)-(|66D and then put the solution obtained into the expression for the energy (57). If we 
suppose that e" > for n > 1 (which can be checked after all by analysing the corresponding 
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equations for the e'^'s) in the sum of the integral equations (p5D-(|6q), only one term will 
contribute in the limit T ^ 0, and we have 



(68) 
where 



e^.«(A) = e^'- + ] dA/i(A-A)el(A), 
ei(A) = -(/i, + B)- j dXMA - A)e-i'^(A), 



4(A) 







^(A), if ei{A)<0 1,, 
if ei(A) > ' ' 



(A) 



(A), if ei''^(A) < 
if ei(A) > 



Hence, in Sutherland's approach we are lead to the concept of two seas: one for real(non- 
string) solutions, another for rapidity A (spin down electrons). 

Suppose that e^''^(A) = at some point X = Q and = at A = Qa, which are the 
Fermi rapidities. Then we can write the equations for the ground state energy as 



L 



2i+l Q 

^ = l-47r^ J dXfi{X)pl'''{X). 



(70) 



k=l 



The equations for the density functions also simplify 

2j+l 2j+l 



/ r-QA roo\ 
Y.Pv'W = hW-Y.\ + UA/i(A-A)a^ 

k=i k=i V-°° -^QaJ 



p'(A) 



(71) 



and 



2i+i 2j+l Q 

Y.M^)=I1 / ^A/i(A-A)pi''=(A) 



k=l k=l ' ^ 

The Fermi boundaries Q and Q\ are determined from the particle number and magnetisa- 
tion equations 



N 

T 

L 




-Qa roo\ 

+ / 1 dAa^(A) 



Qa 



Qa 



dAa^{A) 



(73) 



Due to a theorem by Lieb and Mattis |^6[, in zero magnetic field Sz should be zero, 
which means that spin down electrons constitute half of all electrons M = N/2 and there 
is not string states in the ground state. 

We see that the equations for the ground state are invariant under transformation of 
the index k, which means that we can arbitrarily distribute the colour over the particles. 
They are as gauge modes in the vacuum of gauge theories. 

In the half filled case Qa = and Q = oo (as easy to see in Lai's form of equations Q) 
and following |5|, |9| one can obtain Eq = 1 — 2 log 2. 

It is also possible , following |2(:], to consider the excitations of the model, introduce 



the so called shift functions and write down equations for them. After that, following |37|, 
one can find the S'-matrix for the excitations of the model. It follows as above, that the 
5-matrix of our model is equal to the S-matrix of ordinary t — J model multiplied by unity 
operators in additional colour space. 
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7 Conclusions 



We have constructed a one dimensional model which is based on the t-J model of strongly 
correlated electrons, but which has an additional quantum group symmetry, ensuring the 
degeneration of the states. We use the Bethe Ansatz technique to investigate this model. 
The equations for density fTinetions, written in the thermodynamic limit , demonstrate that 
the additional degrees of freedom of the model behave as gauge modes. The presence of 
these modes, in our opinion, gives rise a possibility to construct a new type of integrable 
models, if one has an interaction between two models of this kind. Also, different topological 
properties, usually appearing in gauge theories, could be present in our model, and would 
be interesting to investigate. 
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